Localization length and impurity dielectric susceptibility in the critical regime 
of the metal-insulator transition in homogeneously doped p-type Ge 
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We have determined the localization length £ and the impurity dielectric susceptibility Ximp 
as a function of Ga acceptor concentrations (A r ) in nominally uncompensated 70 Ge:Ga just below 
the critical concentration (N c ) for the metal-insulator transition. Both £ and Ximp diverge at N c 
according to the functions £ oc (1 — N/N c )~ 1 ' and Ximp <x (N c /N — respectively, with v = 

1.2 ± 0.3 and ( = 2.3 ± 0.6 for 0.99N C < N < N c . Outside of this region (N < 0.99iV c ), the 
values of the exponents drop to v — 0.33 ± 0.03 and £ = 0.62 ± 0.05. The effect of the small 
amount of compensating dopants that are present in our nominally uncompensated samples, may 
be responsible for the change of the critical exponents at N » 0.99iV c . 



PACS numbers: 71.30.+h, 72.80.Cw 

The metal-insulator transition (MIT) in doped semi- 
conductors is a unique quantum phase transition in the 
sense that both disorder and electron-electron interac- 
tion play a key role.EI Important information about the 
MIT is provided by the values of the critical exponents 
for the zero-temperature conductivity, correlation length, 
localization length, and impurity dielectric susceptibility. 
From a theoretical point of view, the correlation length in 
the metallic phase and the localization length in the in- 
sulating phase diverge at the critical point with the same 
exponent v, i.e., they are proportional to \N/N C — 1|~" in 
the critical regime of the MIT. (N is the dopant concen- 
tration and N c is the critical concentration for the MIT.) 
Since direct experimental determination of v is extremely 
difficult, researchers have usually determined, instead of 
v, the value of /i defined by cr(0) oc (N/N c — l)^ ]r where 
cr(0) is the conductivity extrapolated to T — Ofla It is 
also possible to evaluate /i from finite-temperature scal- 
ing of the form cr(JV,IUx T x f(\N/N c - 1|/T») where 
x/y is equivalent to /^Jjd3 Values of v are thenr.obtained 
assuming v — fi for three-dimensional systems. 

In this work we have determined directly the localiza- 
tion length £ and the impurity dielectric susceptibility 
Ximp in neutron-transmutation-doped (NTD), nominally 
uncompensated 70 Ge:Ga just below N c . The applica- 
tion of NTD to isotopically enriched 70 Ge leads to un- 
surpassed doping homogeneity and precisely controlled 
doping concentration. As a result, we have been able 
to approach the transition as close as 0.999iV c from the. 
insulating side and 1.0004A r c from the metallic side.cl 
In zero magnetic field, the low-temperature resistivity 
of the samples is described by variable-range hopping 
(VRH) conduction within the Coulomb gap.B Magnetic 
field and temperature dependence of the resistivity are 
subsequently measured in order to determine directly £ 
and Ximp in the context of the VRH theory B 

This kind of determination of £ and Ximp was per- 



formed for compensated Ge:As by Ionov et alB They 
found £ oc (1 - N/N c )- U and Ximp oc (N c /N - 1)~ ? with 
v = 0.60 ± 0.04 and C = 1.38 ± 0.07, respectively, for 
N < 0.96 N c . The significance of their result is the exper- 
imental verification of the relation 2v w C, that had been 
predicted by scaling theories .llErlid However, the critical 
exponents of compensated samples are known to be dif- 
ferent from those of nominally uncompensated samples. 
Therefore, the present work which probes £ and Ximp 
in nominally uncompensated samples is relevant for the 
fundamental understanding of the MIT. The-previous ef- 
fort to measure Ximp has also contributed.li3~li-3 Hess et 
al. fpund £ = 1.15 ± 0.15 in nominally uncompensated 
Si:P.£3 Since /i w 0.5 was determined for the same series 
of Si:P samples, 2fi w £ was again valid. Katsumoto has 
found ( « 2 and /i«l foe-compensated Al x Gai_ x As:Si, 
i.e., again, 2/jk( applies.Eil Thus, in these cases the con- 
clusion 2v ss £ was reached indirectly, by assuming \i = v. 
The work reported here, on the other hand, determines 
v directly, i.e., we do not have to rely on the assumption 
fx = v in order to study the behavior of £ near N c . 

All of the 70 Ge:Ga samples used in this study were 
prepared by NTD of isotopically enriched 70 Ge single 
crystals. We use the NTD process since it is known to 
produce the most homogeneous dopant distribution.a De- 
tails of the sample preparation and characterization are 
described elsewhere.0 In this study, we determined the 
low-temperature (0.05 — 0.5 K) resistivity of nine sam- 
ples in weak magnetic fields (< 0.4 T) applied in the 
direction perpendicular to the current flow. 

The electrical conduction of doped semiconductors on 
the insulating side of the MIT is often dominated by VRH 
at low temperatures. The temperature dependence of the 
resistivity p(T) for VRH is written in the form of 

P (T)= Po (T)eM(T /TT}, (1) 
where p = 1/2 for the excitation within a parabolic- 
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shaped energy gap (the Coulomb gap) ,u and p = 1/4 for a 
constant density of states around the Fermi level.cJ In our 
earlier work|a we reported that p =1/2 for N < 0.991./V C 
(N c = 1.860 x 10 17 cm -3 ) and that p decreases rapidly as 
N approaches N c from 0.991iV c and becomes even smaller 
than 1/4 when we neglect the temperature variation of 
po{T). However, the variation contributes greatly to the 
temperature dependence of p(T) near N c because the fac- 
tor Tq/T in the exponential terms become very small, 
i.e., the temperature dependencies of po(T) and that of 
the exponential term become comparable. Theoretically, 
Po(T) is expected to vary as po oc T~ r but the value of 
r including the sign has not been derived yet for VRH 
with both p = 1/2 (Rcf. |6l) and p = 1/4 (Ref. [t7|). 
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FIG. 1. Resistivity multiplied by T 1/a vs (a) T~ 1/2 and (b) 
T" 1/4 for 70 Ge:Ga. From top to bottom in units of 10 17 cm" 3 , 
the Ga concentrations are 1.848, 1.850, 1.853, and 1.856. 

Recently, we have shown that the temperature depen- 
dence of the conductivity of the same series of 70 Ge:Ga 
samples within ±0.3% of N c is proportional to T 1/3 at 



0.02 - 1 K.i Since both the T 1/3 dependence of the con- 
ductivity and the Efros-Shklovskii VRH are results of the 
electron-electron interaction in disordered systems, they 
can be expressed, in principle, in a unified form. More- 
over, the electronic transport in barely metallic samples 
and that in barely insulating samples should be essen- 
tially the same at high temperatures so long as the in- 
elastic scattering length and the thermal diffusion length 
are smaller than, or at most comparable to the corre- 
lation length or the localization length. So, the tem- 
perature dependence of conductivity at high tempera- 
tures should be the same on both sides of the transi- 
tion. Such behavior is confirmed experimentally in the 
present system,El i.e., the conductivity of samples very 
close to N c shows a T 1/3 dependence at T« 0.5 K, irre- 
spective of the phase (metal or insulator) to which they 
belong at T = 0. Based on this consideration we fix 
r = 1/3. Figure [j] shows pT r with r = 1/3 for four 
samples (N/N c = 0.993, 0.994, 0.996, and 0.998) as a 
function of (a) T~ 1/2 and (b) T~ 1/4 . All the data points 
lie on straight lines with p = 1/2 in Fig. [j](a) while they 
curve upward with p = 1/4 in Fig. [l](b). This depen- 
dence is maintained even when we change the values of r 
between 1/2 and 1/4. Thus we conclude that the resis- 
tivity of all samples for N up to 0.998-/V c is described by 
the VRH theory where the excitation occurs within the 
Coulomb gap, i.e., Eq. (|l]) with p = 1/2. 

Based on these findings, we evaluate Tq in Eq. (|l|) 
with p = 1/2 and r = 1/3, and show it as a function 
of 1 — N/N c in Fig. |^. The vertical and horizontal er- 
ror bars have been estimated based on the values of To 
obtained with r = 1/2 and r = 1/4, and the values of 
l-iV/(1.858xl0 17 cm- 3 ) and 1- A7(1.861 x 10 17 cm- 3 ), 
where 1.858 x 10 17 cm -3 is the highest concentration in 
the insulating phase and 1.861 x 10 17 cm -3 is the lowest 
in the metallic phase, respectively. 
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FIG. 2. T determined by p(T) oc T _1/3 exp[(T /T) 1/2 ] as 
a function of the dimensionless concentration 1 — N/N c . 
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According to theoryja To in Eq. ([!]) is given by 

k B T « 2.8e 2 /4^e e(JV) £(JV) (2) 

in SI units, where e(N) is the dielectric constant. Here, 
we should note that the condition T < T is needed 
for the theory to be valid, i.e., Tq has to be evaluated 
only from the data obtained at temperatures low enough 
to satisfy the condition. This requirement is fulfilled 
in Fig. H for all the samples except for the one with 
N = 0.998A C . Concerning this latter sample, we will 
include it for the determination of £ and Ximp (Fig. |^) 
but not for the calculation of the critical exponents. 




FIG. 3. Logarithm of resistivity vs B 2 at constant temper- 
atures for the sample having N = 1.840 x 10 17 cm -3 . From 
top to bottom the temperatures are 0.095 K, 0.135 K, and 
0.215 K, respectively. The solid lines represent the best fits. 

Our next step is to separate To into e and £. For 
£/A <C 1, the magnetoresistance is expressed as 

\n[p(B, T)/p(0, T)] a 0.0015 (£/A) 4 (T /T) 3 / 2 , (3) 

where A = \[K[eB is the magnetic length in SI units.i 
According to Eq. (||) , the magnetic-field variation of In p 
at T = const, is proportional to B 2 , i.e., \np(B,T) — 
lnp(0, T) + C(T)B 2 , and the slope C(T) is proportional 
to J 1-3 / 2 . In order to demonstrate that these relations 
hold for our samples, we show for the N = 0.989iV c sam- 
ple In p(B, T) vs B 2 in Fig. | and C(T) determined by 
least-square fitting of d In p/dB 2 vs T" 3 / 2 in Fig. [l| Since 
Eq. (|J) is equivalent to 

7 = C(T)/T- 3 / 2 « 0.0015 (e/?i) 2 £ 4 T Q 3/2 , (4) 
£ is given by 

e«5.1(Ve) 1/2 7 1/4 T - 3/8 . (5) 

In this way we have determined 7 for nine samples. The 
inset of Fig. shows 7 as a function of Tq. The value of 7 



is almost independent of To, and if one assumes 7 oc T Q S , 
one obtains a small value of 5 = 0.094 ± 0.005 from least- 
square fitting. Figure J| shows £ and Ximp = e — e/i 
determined from Eqs. (g) and (|5j). Here, e/j is the di- 
electric constant of the host Ge, and hence, Ximp is the 
dielectric susceptibility of the Ga acceptors. We should 
note that both £ and Ximp are sufficiently larger than the 
Bohr radius (8 nm for Ge) and eh = 15.4 (Ref. |l2|), re- 
spectively. According to the theories of the MIT, both 
£ and Ximp diverge at iV c as £(iV) oc (1 — N/N c )~" and 
Ximp(-^V) K {N c /N— 1) _? , respectively. We find, however, 
both £ and Ximp do not show such simple dependencies 
on N in the range shown in Fig. |^, and that there is a 
sharp change of both dependencies at N « 0.99A C . On 
both sides of the change in slope, the concentration de- 
pendence of £ and Ximp are expressed well by the scaling 
formula as shown in Fig. |^. Theoretically, the quanti- 
ties should show the critical behavior when N is very 
close to N c . So f = 1.2 ± 0.3 and £ = 2.3 ± 0.6 may 
be concluded from the data in 0.99 < N/N c . However, 
the other region (0.9 < N/N c < 0.99), where we obtain 
v = 0.33 ± 0.03 and £ = 0.62 ± 0.05, is also very close to 
N c in a conventional experimental sense. 
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FIG. 4. Slope dlnp/dB 2 vs T" 3/2 for the sample having 
N = 1.840 x 10 17 cm" 3 . The solid line represents the best fit. 



As a possible origin for the change in slope, we refer 
to the effect of compensation. Although our samples are 
nominally uncompensated, doping compensation of less 
than 0.1% may present due to residual isotopes that be- 
come n-type impurities after NTD. In addition to the 
doping compensation, the effect known as "self compen- 
sation" may play an important role near N C S3 It is em- 
pirically known that the doping compensation affects the 
value of the critical exponents. Rentzsch et al. stud- 
ied VRH conduction of n-type Ge in the concentration 
range of 0.2 < N/N c < 0.91, and showed that T van- 
ishes as T oc (1 - N/N c ) a with a « 3 for K = 38% 
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and 54%, where K is the compensation ratio.0 Since 
a w v + ( [Eq. (§)], we find for our 70 Ge:Ga samples 
a = 3.5 ± 0.8 for 0.99 < A/A c < 1 and a = 0.95 ± 0.08 
for 0.9 < A/A c < 0.99. Interestingly, a = 3.5±0.8 agrees 
with a m 3 found for compensated samples. Moreover, 
we have recently proposed the possibility that the con- 
ductivity critical exponent /i « 1 in the same 70 Ge:Ga 
only within the very vicinity of N c (up to about +0.1% 
of N C )E An exponent of fi = 0.50 ± 0.04, on the other 
hand, holds for a wide region of N up to 1.4iV c .tj Again, 
fi w 1 near A c may be viewed as the effect of compen- 
sation. Therefore, it may be possible that the region 
of N around N c where v « 1 and fi » 1 changes its 
width as a function of the doping compensation. In the 
limit of zero compensation, the part which is charac- 
terized by v k. 1 and fi ss 1 vanishes, i.e., we propose 
= 0.33 ± 0.03, C = 0.62 ± 0.05, and fi = 0.50 ± 0.04 for 
truly uncompensated systems and that Wegner's scaling 
law of v = fi is not satisfied. In compensated systems, 
on the other hand, Wegner's law may hold as it does in 
the very vicinity of N c . The experiment on compensated 
ALjGai-zAsiSi that showed ( « 2 and fi « 1 (Ref. [hi]) is 
also consistent with the law. However, our preceding dis- 
cussion needs to be proven in the future by experiments 
in samples with precisely and systematically controlled 
compensation ratios. 
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FIG. 5. Localization length £ vs 1 — N/N c (lower data set) 
and the dielectric susceptibility Ximp arising from the impu- 
rities vs Nc/N — 1 (upper data set). The inset shows the 
coefficient 7 defined by Eq. (Q) as a function of To- 
ll! summary, we have determined directly the local- 
ization length and the dielectric susceptibility arising 
from the impurities in nominally uncompensated NTD 
70 Ge:Ga samples near the critical point for the MIT. 



While the relation 2i/k( predicted by scaling theorjJilrtLil 
holds for 0.9 < N/N c < 1, the critical exponents for lo- 
calization length and impurity susceptibility change at 
N/N c k 0.99. The small amount of doping compensa- 
tion that is unavoidably present in our samples may be 
responsible for such a change in the exponents. 
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